Giant Nernst-Ettingshausen Oscillations in Semiclassically Strong Magnetic Fields 
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We consider the Nernst-Ettingshausen (NE) effect in the presence of semiclassically strong mag- 
netic fields for a quasi-two-dimensional system with a parabolic or linear dispersion of carriers. We 
show that the occurring giant oscillations of the NE coefficient are coherent with the recent exper- 
imental observation in graphene, graphite and bismuth. In the 2D case we find the exact shape 
of these oscillations and show that their magnitude decreases/increases with enhancement of the 
Fermi energy for Dirac fermions/normal carriers. With a crossover to 3D spectrum the phase of 
oscillations shifts, their amplitude decreases and the peaks become asymmetric. 

PACS numbers: 72.15.Jf, 72.20.Pa 
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The Nernst-Ettingshausen (NE) effect in metals [l| is a 
thermoelectric counterpart of the Hall effect. The effect 
consists in induction of an electric field E y normal to the 
mutually perpendicular magnetic field H (|| z) and tem- 
perature gradient S7 X T. All electric circuits are supposed 
to be broken: J x = J y = and heat flow along y-axis 
to be absent (adiabatic conditions). Quantitatively, the 
effect is characterized by the NE coefficient. 
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The NE coefficient varies by several orders of magnitude 
in different materials ranging from about 7mV ■ K~ 1 T~ 1 
in bismuth up to W~ 5 mV • K~ 1 T~ 1 in some metals Q. 

The NE effect was discovered in 1886 and remained 
poorly understood until 1948 when Sondheimer us- 
ing the classical Mott formula for the thermoconductivity 
tensor, calculated v for a degenerated electron system. It 
has been linked to the energy derivative of the Hall an- 
gle 6 = a xy /(7 xx . Within this model, v was found to eter 7 |12|. Il 
be independent on the magnetic field in weak fields and 
to decrease as H~ 2 in the region of semiclassically strong 
fields, where the cyclotron frequency uj c is larger than the 
inverse scattering time r" 1 . In 1964, Obraztsov Q sug- 
gested that magnetization currents (i.e. electric currents 
induced due to inhomogeneous distribution of magneti- 
zation in the sample) can contribute supplementary to 
the NE effect. 

The giant oscillations of v were firstly experimentally 
observed in 1959 in zinc by Bergeron et al [5j who quali- 
tatively ascribed the phenomenon to crossing of the elec- 
tronic Fermi energy by Landau levels (LL). Similarly to 
de Haas - van Alphen (dHvA) oscillations of magnetiza- 
tion and Shubnikov - de Haas (SdH) oscillations of con- 
ductivity, in the NE oscillations the corresponding quan- 
tizing fields are given by Lifshitz-Onsager condition 0] : 



where S (/x) is the cross section of Fermi surface (FS) of 
the orbital electron motion at p z = 0, fi is the chemical 
potential, k is integer. Here 7^ = 7+57^-c with a = ±1, 
and the electron cyclotron mass m* = @. 

Very recently, the NE effect has been measured 0, H| 
and theoretically analyzed 0] in graphene. Surprisingly, 
it has been found that v changes its sign at quantizing 
field in graphene while it has maxima in zinc [5[ and bis- 
muth [T(| ■ Zhu et al. [ll| demonstrated that such untyp- 
ical behavior of v(H ) observed in graphene is not repro- 
duced in graphite. They concluded that piling of multi- 
ple graphene layers leads to a topological phase transition 
in the spectrum of charge carriers, so that graphite be- 
haves as a 3D crystal despite of its apparent structural 
anisotropy and of similarity of its electronic properties 
to those of graphene. 

Another challenging property of quantum oscillations 
is the possibility to distinguish between two types of 
charge carriers, having the topologically different param- 
for the normal carriers (NC) with 
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parabolic 2D dispersion and linear LL quantization: 



NC: e(p ± ) 



2m 1 



m ( 1 
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and 7 = for the Dirac fermions (DF) having the linear 
two-branch spectrum and ~ fc 1 / 2 LL quantization: 



DF: 



±v\p±\, e k = ± [4mv 2 n B H k] 
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(1) 



p± and m± being momentum and effective mass in the 
plane normal to the magnetic field, m is the free electron 
mass, v is the Fermi velocity and \i B = eh/2mc is the 
Bohr magneton. 

In this Letter we propose a simple thermodynamic 
approach to the description of the NE effect which al- 
lows linking the oscillations of the NE coefficient to the 
oscillations of the magnetization. Both thermal (Sond- 
heimer) and magnetization (Obraztsov) contributions to 
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the Nernst coefficient are evaluated analytically for a 
quasi-two dimensional (q2D) electronic system with ei- 
ther parabolic or Dirac spectrum. In the 2D limit for the 
Dirac spectrum we recover the behavior of the NE coef- 
ficient observed in graphene @, [1] while the recent data 
of Zhu et al. on graphite correspond to the 3D limit. 

Thermodynamic approach. The NE coefficient is mea- 
sured in the absence of the electric current flowing 
through the system along the temperature gradient. This 
is why the system can be assumed to be in thermody- 
namic equilibrium where the electrochemical potential 
p + e<p = const, with (p being the electrostatic potential. 
Hence the effect of the temperature gradient is reduced 
to the appearance of an effective electric field in the x- 
direction E x = \7 x /j/e. In this way, the problem is re- 
duced to the classical Hall problem, which allows us to 
obtain the thermal contribution to the NE coefficient: 
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where a xx is the diagonal component of the conductivity 
tensor, n is the concentration of carriers. This simple for- 
mula reproduces Sondheimer's result for a normal metal, 
fluctuation contribution to the NE coefficient in a super- 
conductor above T c , etc. [3, 15 1. 

The additional contribution to the NE coefficient ap- 
pearing due to the spatial dependence of magnetization 
in the sample can be found from the Ampere law. The 
magnetization current density is j ma s = x B, where 
B = H + 47rM, H is the spatially homogeneous exter- 
nal magnetic field, M is the magnetization, which can be 
temperature and, henceforth, coordinate dependent. In 
the case under consideration one can express the mag- 
netization current as j™ ag = -c(dM/dT)V x T [| and 
the corresponding contribution to the electric field in the 
y- direction (Nernst field) as £™ ag = p yy jy 



mag wncrc 



p yy is the diagonal component of the resistivity tensor 
(■Pyy = Pxx)- The magnetization contribution to the NE 
coefficient reads as 



,mag _ CPyy ( dAl\ 

H \dT ) 



(3) 



evaluate these quantities no supplementary knowledge of 
the transport properties of the system is needed. These 
derivatives can be expressed in terms of the thermody- 
namic potential of the system: 
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(4) 



To be more specific, we consider the quasi-2D system 
with the dispersion 



2t sin ^-d. 



(5) 



This model allows us to describe the 2D-3D dimensional 
crossover by variation of the hopping parameter t from 
t-2D = to ~ £f- The corresponding expression for 
the oscillating part of Q (denoted by tilde), derived by 
Champel and Mineev for the parabolic dispersion [16[ 
(see also [l?}) and generalized in [l8| for the arbitrary 
e±(p±) reads: 
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with ip(Xl) = 
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l=l,a=±l 



and 
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Re$, CT (p, H) , (6) 



<5> la (p, H) = Jo[ 2ttZ 
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Here fcs = 1, A = 2j-f-, T is the Dingle LL broadening 
and Jo is the Bessel function. We present Eq. ((6]) in the 
most general form using the parameters S(p) at p z = 0, 
m*,w c and j a . For NC S = 27rm^/i, m* = m±, uj c = 



±Tc aild y = 2 • 2 



1 m 



2^cr; for DF 5 = tt^, 



TO = fr, 
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a. In the present derivation 
we assume a Lorentzian broadening of Landau levels with 
a constant T. Such approximation can be justified for 
w c <C £j? in the case of 3D system. In 2D systems it is 
expected to be valid only in the low field regime w c < 
t _1 . The oscillating parts of the chemical potential and 
magnetization can be expressed using Eq. as: 



The Eqs. © and ^ reveal the essential physics of 
Nernst oscillations in the quantizing magnetic fields. In 
particular, one can see that the NE coefficient is depen- 
dent on the diagonal components of conductivity and 
resistivity tensors. Their oscillations as a function of 
the magnetic field constitute the SdH effect. The giant 
Nernst oscillations have been observed even in the regime 
where the SdH effect is weak in graphene (at H < 3T) 



and in graphite [11 



This is why one should attribute 
the giant NE coefficient oscillations to the remaining fac- 
tors in the Eqs. ([2]) and ((3]), namely, to the tempera- 
ture derivatives of the chemical potential and magneti- 
zation, dp/dT and dM/dT, respectively. Remarkably, to 



dp ImS' 1 ' dM n dp 

dT = ~l + 2ReSM' ~dT ~ HdT' 
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(8) 
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Z = 1,(T = ±1 



and (x) is the derivative of the order of a = 0, 1 of 
the function ifi. One can see from Eqs. ([3]) and (JSJ) that 
the NE coefficient oscillates proportionally to the deriva- 
tive of magnetization over temperature. This shows an 
important link between NE and dHvA oscillations, which 
is universal and independent on the dimensionality of the 
system and of the type of carriers. 
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It is convenient to express the NE coefficient as 

v (H) + v(H) (10) 



v = v theim + ^ mag 



with vq (H) and v (H) being the background and oscil- 
lating parts. The background part can be evaluated in 
the Drude approximation as [l5| 



u (H) = 



7T 2 T (T 



6m*c \sfJ 1 + (w c r) 



(11) 



The account for magnetization currents leads to the cor- 
rection of the order of (eft)~ 2 with respect to Sond- 
heimer result described by Eq. (fTTj) . 

The oscillating part of the Nernst coefficient can be 
written using Eqs. I©,© and (© as: 



u(H) = -2ttk (H) 



ImSi 1 } 
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with 



K (H) 
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e 2 nc H 2 



(12) 



(13) 



In the Drude approximation for NC 

r 11 + 2(w c t) 5 



K Drudc 



m*c (lu c t) 2 1 + (w c t)' 



(14) 



Equation (|T2"j) describes oscillations of the NE effect in 
the most general form. It is valid for any type of the 
dispersion e±(p±) if T, t •< /i. 

The 2D case: graphene. We start analysis of the Eq. 
(fl"2"|) from the pure 2D case where t = 0. In the low- 
temperature limit 2ir 2 T < hu c in Eq. jB) A < 1, hence 
ip(Xl) 1 - \\ 2 l 2 . For m* < 0.02m and H = 10T 
(typical in graphene experiments) this yields T < 10K. 
Since m* <C m we neglect also the Zeeman splitting, 
assuming that r y a = 7 = for NC and r ) a = 7 = 5 for 
DF. The series S^ ^ and S^ 1 ! in Eq. (|12[) in this case can 
be summed exactly which gives: 
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In the experimental configuration corresponding to the 
measurement of the NE effect in graphene, the number 
of particles n is fixed, so that [la ]: 



V 5 M / HjT 



(2nh) 2 V ^ )h: 



const 



(16) 



(we assume the volume V = 1). This relation implicitly 
determines the dependence of /1 on H, T for the given n. 
We note that the chemical potential /1 itself is a function 
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FIG. 1: Normalized Nernst-Ettingshausen (NE) oscillation as 
function of the inverse magnetic field and carriers concentra- 
tion for normal carriers (NC) and Dirac fermions (DF). De- 
pendence v{H~ ) for DF has the same profile as for NC but 
shifted on half period. Vertical lines shows the quantization 
condition |T}. 



of H as follows from Eq. (|16[) . which in the 2D case can 
be written as: 



S{n) 



h 2 



■ arctan ■ 



(2tt/i) 

This equation can be inverted for S(u) 
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Equation (|T7l) yields the dependence /i(n,H). Substitut- 
ing it to Eq. (|15l) after some cumbersome algebra one 
can find the oscillating part of the Nernst coefficient ex- 
plicitly: 



(n,H) 



2ir 3 T k(H) . „ / hen 

tt-^ft srnz7r n 7 

3 hw c sinh V e H 



huj c 



(18) 

that is a strongly oscillating function. It crosses zero 
at the intersections of LL and chemical potential, given 
by the condition H = Hk a defined by (fl]). The field 
depended factor k (H) is governed by magnetoresistance 
and is given by Eq. tfTB")) . At lj c t < 1 where SdH oscilla- 
tions are small, k (H) can be roughly estimated using the 
Drude approximation (|14|) . In particular, approaching 
the limit w c r ~ 1 and assuming T ~ h/2r we obtain that 
k (H) ~ ^j— and the amplitude of NE oscillations is gi- 
ant in comparison with the background: i^ 2D ^ ~ ^f-^o- 
At higher fields uj c t > 1, in the quantum Hall regime, 
the shape of oscillations of the NE coefficient is affected 
by strong variation of the magnetoresistance and Dingle 
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temperature. This can be taken into account by substitu- 
tion of the field dependent magnetoresistance and Dingle 
temperature into Eqs. (|T2" ]) .(fT3" ]) . 

The given by Eq. ([TH]) profiles of 2D NE oscillation 
as function of H and n for DF and NC are presented 
in Fig.l. Both our theory for DF and experiment in 
graphene @, Q show a sin-like profile of the signal whose 
amplitude slightly decreases with increasing n. This ten- 
dency contradicts to the earlier theoretical predictions of 
the classical Mott formula Q that has been derived for 
a Boltzmann gas of electrons. In contrast, the amplitude 
of NE oscillations increases with increasing n for the NC 
in a qualitative agreement with the Mott formula. 

Quasi-2D and 3D cases. In order to describe the NE 
effect in the general quasi-2D case where t ^ the Bcsscl 
function in the Eq. jTJ should be taken into account. 
The sums (J9j) can be reduced to the integrals by means 
of the Poisson transformation. Then integration can be 
done analytically resulting in 
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(19) 



(20) 



where 6 ka (H) = £- - (H^ - H ka , 



" n The NE 

coefficient is obtained by substitution of the Eqs. (fl"9|) . 



and (J20J) to Eq. ([12]). Resonances at iS k a (H) = ±-£- 
in v (H) appear when the chemical potential crosses the 
quantized slices of maximal (minimal) cross sections of 



the corrugated cylinder FS S, 



max(min) 



S ± 47TtTO* 



In the wide quasi 2D interval t < (ftw c ) /Y the be- 
havior of v^ q2D ^ (iJ) close to H = Hka can be studied 
selecting in (fT9|) and (|20| only the resonant terms. With 
growth of t the positions of zeros shift from Im5k& (H) = 
to ImSka {H) = ±t^- The superposition of two (for 
S'max and ^min) series of resonances leads to the beats in 
v (H) oscillations. 

In the 3D limit t > (huj c j 2 /V , RcS {0} < 1, so that 
can be neglected in the denominator of Eq. (fl"2"|) . In 
the vicinity of H = H^ a one finds 



(H) = T 



tx Tk(H) 



12 



(tfuv c ) 



1/2 



Re- 



3/2 ' 



(21) 



We assumed here the constant /z and neglected Zeeman 
splitting, taking 5k.±i = <5fc- The resonances in v (H) 
described by Eq. (f2~Tj) have the form of asymmetric spikes 



with 
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3.4 as shown in FigHJ In 



the Drudc approximation, the amplitude 
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v Q (H) 



(22) 



is giant if > 1. 

For 2D systems our calculations arc valid for mag- 
netic fields uj c < t -1 where one can neglect the quan- 
tum Hall oscillations of conductivity. At higher fields 
the approach of Girvin and Jonson [lj|, based on the 
generalized Mott formula for the thermopower tensor for 
2D systems, seems to be more relevant. In 3D case the 
range of applicability of our theory is given by u c < £p. 
Recently Bergman and Oganesyan Q extended the ap- 
proach of Ref. [l9[ to calculate the off-diagonal thermo- 
electric conductivity a xy for a 3D system at uj c ~ Ef- 
Although a xy constitute only the part of the NE coeffi- 
cient v = — (p xx Oi xy + p X yCtyy) JH , they reproduce quite 
well the measured in graphite jll| sawtooth dependence 
of v(H), having the characteristic {Hk — H)~? divergen- 
cies at resonances. 

In conclusion, we have obtained an analytical expres- 
sion for the oscillating NE constant in a 2D system with 
an arbitrary electron dispersion, describing the recent ex- 
perimental results in graphene and predicting a qualita- 
tive difference in the NE oscillations for NC and DF. We 
show that the giant oscillations of the NE coefficient pre- 
dicted and observed in a 2D case (graphene) decrease 
significantly as the spectrum acquires a 3D character 
(graphite). We describe analytically the shape of NE 
oscillations. The NE oscillations are proportional to the 
temperature derivative of the dHvA oscillations. 
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